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Abstract 



We discuss two-photon and hadronic production of Be mesons in nonrelativistic bound state approx- 
imation and to lowest order in the couphng constants a. and a^. It is shown that in photon-photon 
cohisions, heavy quark fragmentation is dominated by recombination of b and c quarks up to the 
highest accessible transverse momenta. In contrast, in hadroproduction, which at high energies 
mainly involves gluon-gluon collisions, the fragmentation mechanism dominates at transverse mo- 
menta pt > "mBci providing a simple and satisfactory approximation of the complete 0(af) results 
in the high-pT regime. Contradictions in previous publications on hadroproduction of mesons are 
clarified. We also present predictions for cross sections and differential distributions at present and 
future accelerators. 



^Talk presented by K. Kolodziej at The Second German-Polish Symposium "New Ideas in the Theory of 
Fundamental Interactions", Zakopane, Poland, 11-15 September, 1995. 

^Supported by the German Federal Ministry for Research and Technology BMBF under contract 
05 6MU93P, and by the EC-program HCM under contract CHRX-CT93-0132 

^Supported by the Polish Ministry of Education under contract No. PB 659/P03/95/08 



Since the top quark is too short hved for the formation of quarkonium-hke resonances, Be 
mesons are most probably the only flavoured heavy quark resonances in nature. Because of 
flavour conservation in strong and electromagnetic interactions, the Be ground state must 
decay weakly. The nonrelativistic nature of these bound states provides unique possibilities to 
compute genuine nonperturbative quantities such as fragmentation functions and weak matrix 
elements, and to study interesting aspects of the strong and weak dynamics of hadrons. 

A limit on Be production has been reported recently by the CDF-collaboration at the Tevatron 
[Q]. At the LHC, the production rates are predicted to be large enough for a detailed study 
of the production and decay properties [^. Also at linear colliders in the TeV energy range. 
Be mesons produced in collisions of Compton or bremsstrahlung photons may come into 
experimental reach |^. 

In this talk, we report about two recent studies [|, Q of Be production. We discuss the frag- 
mentation and recombination mechanisms and compare the relative importance of them in 
photon-photon and hadron-hadron collisions. As a main result, we clarify quantitatively the 
validity of the hard scattering description in terms of heavy quark fragmentation functions. 
Furthermore, on the basis of two completely independent calculations, we resolve contra- 
dictions in previous publications on hadronic production of Be mesons Finally, we 
present the most relevant integrated cross sections and differential distributions as predicted 
in lowest -order perturbation theory and nonrelativistic bound state approximation. 

In photon-photon collisions. Be mesons are produced in association with b- and c-quark jets: 

77 Bebc. (1) 

This is also the case in hadronic collisions, where gluon-gluon scattering, 

99 ^ Bebc, (2) 

is the dominant subprocess at high energies. In general, one can distinguish two production 
mechanisms, namely heavy quark fragmentation and recombination. We have found these 
mechanisms to contribute quite differently in the two reactions, which makes a comparison 
very interesting. 

The twenty 0(a^Q;^) Feynman diagrams of process (|^) can be classified in three gauge invariant 
subsets characterized in Fig. 1: subset (I;,) with the 6-quark line coupled to the primary 
photons, subset (Ic) with the c-quark line coupled to the primary photons, and subset (II) 
with the 6-quark line coupled to one of the primary photons and the c-quark line to the other 
one. The three subsets can be interpreted physically as describing, respectively, bb production 
and subsequent b fragmentation, b BeC, cc production and subsequent c fragmentation, 
c — > Beb, and simultaneous production of a bb and cc pair and recombination of the b- and 
c-quark into a Be meson. The gluon-fusion process (H) involves a larger number of Feynman 
diagrams because of the presence of the gluon self-coupling. To one has thirty-six 

diagrams in total. More importantly, they cannot be divided up in gauge invariant subsets 
corresponding to the production mechanisms typified in Fig. 1. 

In the following, we concentrate on the production of the pseudoscalar (Be) and vector (B*) 
bound states. Furthermore, we consider the nonrelativistic limit, in which the relative momen- 
tum of the constituents and their binding energy are neglected relative to the b- and c-quark 
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masses. In this limit, the masses of the bound states are equal to the sum of mb and rric, 
and the momenta of both constituent quarks are proportional to the bound state momentum. 
Moreover, the amplitudes for the production of S- waves factorize into hard scattering ampli- 
tudes for 77 (or gg) — > bbcc times the S'-wave function at the origin. The latter can be related 
to the B^*^ decay constants. All this results in the simple substitution rule P, ^ 

vimMPc) = ^(]/ - M)Uss, , (3) 

indicated in Fig. 1 by the black blob, where v{pi) and u{pc) denote the b- and c- quark spinors, 
respectively, and Ussz = l^if) is the spin projector for Be (B*). Note that the colour structure 
is not accounted for in Eq. (0). 




Fig. 1. Characteristic topologies of the lowest- 
order Feynman diagrams contributing to 77 
BM. 



Fig. 2. Transverse momentum distributions of 
Be mesons produced in 77 Bcbc at ^/s = 100 
Ge V for the different production mechanisms 
and approximations described in the text. 



The squared matrix elements for the processes (0) and (H) have been calculated independently 
using two different methods, the traditional trace technique and the method of helicity ampli- 
tudes. The results are found to be in perfect agreement. In addition, we have checked external 
gauge invariance, that is the vanishing of the matrix elements when the polarization vector of 
any of the initial photons or gluons is substituted by its momentum. For process (|I|), we have 
also tested internal gauge invariance, that is the independence of the matrix element of the 
gauge parameter in the gluon propagator. Finally, we have double-checked the phase space 
integration using two different Monte Carlo routines. For more details of the calculation we 
refer to refs. and 

In addition to the above Oi^a^al) and calculations, we have also studied the factor- 

ized description of the processes (|1]) and (0) in terms of bb{cc) production followed by 6(c) 
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fragmentation |]T2 



das, = da^i ® Di^^A^) i^^^cc ® D^^bA^))- (4) 

The relevant fragmentation functions Di{z) and D^{z) have been derived from perturbation 
theory They are known to provide a perfect approximation of the energy distribution 
dcx/dz in e^e^ Bccb |TD[ in order a^a^. Clearly, in photon-photon and hadron-hadron 
production such a factorized description cannot be expected to work close to threshold and 
at small pt, where the quark masses play a role. The question is how well this approximation 
works at high-p^^. Previous studies of this issue only give qualitative and partly contradicting 
answers This motivated us to compare the factorized approximation (|^) with our com- 

plete lowest -order calculations, and to determine the region of validity of (||) quantitatively. 

The numerical results plotted in Figs. 2-6 have been obtained with the following values of the 
parameters: 

mb = 4.8GeV, m, = 1.5 GeV, a = 1/129, /^^ = J^, = 0.4 GeV [11]. (5) 

Furthermore, we have used the running coupling constant as{Q^) in leading logarithmic ap- 
proximation for five flavours and normalized to as(m|) = 0.113. Additional specifications are 
given when needed. 

Fig. 2 illustrates transverse momentum distributions of the B^. in 77-production (|1|). We see 
that the description in terms of bb production and 6-fragmentation and the corresponding 
subset (If,) of the diagrams of Fig. 1 give distributions (dotted and dashed curves, respectively, 
labeled by b) which are very similar in shape and normalization except in the low-p^ region, 
where the factorized approximation is expected to break down. As far as the shape is concerned 
this is also true for cc production and c fragmentation (curves labeled by c). However, in 
this case the approximation (H) fails to reproduce the correct magnitude of the cross section 
predicted by the diagrams subset (Ic) of Fig. 1. Note that primary c-quark production is 
enhanced to 6-quark production by a factor 16 due to the ratio {Qc/ QbY of the electric charges. 
On the other hand, the radiation of a cc-pair from a 6-quark leads to a harder pt spectrum for 
the Bf. bound states than the radiation of a 66-pair from c-quarks. Finally, the most important 
observation is that the recombination mechanism, represented by the diagrams (II) of Fig. 1, 
dominates B^. production not only at low pt, as one could have expected, but also in the 
high-pT region up to the kinematical limit. In other words, the familiar description of high-pr 
hadron production in terms of the production and fragmentation of quarks is inadequate for 
single Be production in 77-scattering. 

The main features of the transverse momentum distributions of the B^. produced in the sub- 
process gg Bcbc are illustrated in Fig. 3. Here, we see that the fragmentation description 
(H) indeed approaches the pT^distributions resulting from the complete 0{ag) calculation, but 
only in the tails of the distributions. In order to demonstrate the effect of the quark masses 
on the fragmentation kinematics, we have assumed three different relations between daughter 
and parent momenta: 

Pt = sin Ob, pL = PTCoidb, fi = M{l), m,, (II), (HI), (6) 

where \/3 is the gluon-gluon cm. energy. Case I obeys the physical phase space boundaries. 
The choices II and HI have been considered in refs. [|12| and |T^, respectively. Fig. 3 shows 
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that the mass ambiguities of the fragmentation approach increase as s decreases and that 
they become non-neghgible at < 40 GeV. Comparing the approximations with the 0{aj) 
resuhs, one observes a shght preference for choice I or II. Most interesting, however, is the 
difference to 77 —>■ Bcbc where heavy quark fragmentation is completely subdominant even at 
large pt- Apparently, the presence of the gluon self-coupling and colour factors has a drastic 
influence on the relative importance of the fragmentation and recombination. 
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Fig. 3. Transverse momentum distributions of 
Be mesons produced in gg —i- Bcbc at \/I = 40 
and 100 GeV: complete 0{ag) calculation (cir- 
cles) and approximation eq. (4) (solid curves). 
The labels I-III refer to the kinematics specified 
in eq. (6). 



Fig. 4. Transverse momentum distribution 
of the Be in pp{pp) collisions at the Tevatron 
(LHC) energies: complete 0{ag) calculation 
( circles ) and fragmentation approximation eq. 
(4) (solid curves). The labels I and III refer to 
the kinematics eq. (6). Results are shown with- 
out a rapidity cut (full circles) and for \y\ < 0.5 
(empty circles). 



Predictions for the pt distributions of Be mesons in pp and pp collisions at Tevatron and LHC 
energies are obtained by convoluting the (^f^-subprocess distributions with the MRS (A') gluon 
structure functions |T^. The results with and without cuts in rapidity are plotted in Fig. 4 and 
compared with the fragmentation approximation. We have assumed the scale = p^ + M"^ in 
both as(Q^) and the structure functions. The evolution effects in the fragmentation function 



Di{z) are ignored for consistency. These effects are studied in ref. [0. We see that after 
convolution the 0{aj) calculation and the fragmentation description (^) are in reasonable 
agreement at pr > 10 GeV. This can be understood from the properties of the unfolded pt- 
distributions illustrated in Fig. 3 and from the rise of the gluon density at small x which 
favours contributions from the smallest possible subenergies s, and hence from the tails of 
the spectra. With decreasing pt, the fragmentation picture gradually breaks down and at 
Pt < 5 GeV only the complete 0{a'^) calculation makes sense. Furthermore, the sensitivity to 
the kinematical prescription (|^) decreases slowly with increasing pt- 
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The total hadronic cross sections for Be and B* production are shown in Fig. 5. The typical rise 
of a with energy is due to the rise of the gluon density as x approaches Xmin and the peaking 
of a near threshold. In order to demonstrate the scale dependence of a, we have indicated the 
Be cross sections at ^/s = 0.1, 1 and 10 TeV for different choices of Q^. As one can see, at 
lower energies, the notorious scale ambiguity of leading logarithmic approximations leads to an 
uncertainty of more than one order of magnitude. Only at very high energies, the uncertainty 
shrinks to a factor two. It is interesting to note that the differences in the predictions resulting 
from different parametrizations of the gluon density, e.g. CTEQ2 |T5[ instead of MRS(A') 
would be invisible in Fig. 5. Other uncertainties, connected with the decay constants / 
and the effective quark masses nib and rric, amount at least to another factor of two. 
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Fig. 5. Total cross sections for pp{pp) — > 
versus the cm. energy. The lower 
(upper) curve show the predictions for Be (B*) 
and the scale Q"^ = pj^ + M^. The symbols in- 
dicate expectations for other choices of . 
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Fig. 6. Total cross sections for 77 — > Be be 

versus fixed 77 c. m. energy, and after convolu- 
tion with Compton and brems strahlung spectra, 
versus the e^e~ cm. energy. The lower (up- 
per) curves correspond to Be (B*) production. 
The scale of as is chosen to be p^ + M^. 



The second main motivation for our own studies, besides examining the validity of the factor- 
ized fragmentation approach, has been the wish to resolve the confusion created by previous 
calculations ||^ - |^ of hadroproduction of Be mesons which have contradicted each other. 
Because of the ambiguities and uncertainties pointed out above, comparison was not always 
easy. We have been very careful in adjusting parameters and the gluon density to the choice in 
the respective calculation we considered. We also focused as much as possible on the compar- 
ison of gg cross sections and distributions being most transparent. In cases where gg results 
were not given, we compared the convoluted cross sections. The outcome of this comparison 
was surprising. As described in detail in ref. 0, we were not able to reproduce any of the 
previous numerical predictions. In the meantime, the authors of ref. Q revised their work 
1^ and now agree with our results. Also the authors of corrected a normalization error in 
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(this paper appeared almost simultaneously with M), so that the gg cross section coincide 



with what we found except at very high energies. However, the pt distributions given in 
still disagree substantially with ours. 

In order to evaluate the observability of mesons at the Tevatron and LHC, it is useful to 
integrate the p^^distributions of Fig. 4 over pt > p™". Assuming an integrated luminosity of 
100 pb~^, one can expect about 10^ Be mesons with p™" = 10 GeV at the Tevatron, without 
taking into account contributions from the production and decay of B* mesons and heavier 
states. This rate should be sufficient for discovery in the channel Be J/ipX, for which a 
branching ratio of the order of 10% is predicted Jl^. In fact, first results of such a search have 



already been reported in |1|]. Finally, at the LHC for lOOfb"^, one can expect 10^ direct Be 
mesons at p™*" = 20 GeV. This rate should then allow a more detailed study of production 
and decay properties. 

Predictions on the production rates of Be mesons in photon-photon collisions at future e^e~ 
machines are obtained by folding the total cross sections for 77 — > B^*^bc with the photon 
spectrum generated by Compton back-scattering of high intensity laser light on beams fl^ , 



or with the Weizsacker-Williams bremsstrahlung spectrum |T9[. The expectations are illus- 
trated in Fig. 6, where we have plotted the convoluted cross sections for Be and B* production 
versus the e^e~ centre-of-mass energy, together with the unfolded cross sections as a function 
of fixed 77 centre-of-mass energy. Because of the long soft tail of the Compton spectrum, and 
the shape of the 77 cross section which peaks just above threshold, the convolution increases 
the cross sections substantially for energies above 100 GeV. At a 500 GeV linear collider and 
for an integrated luminosity of 10 fb~^, one can produce about 100 Be and 400 B*. The yield 
of Be mesons from bremsstrahlung photons is invisibly small at LEP energies, but increases 
logarithmically with energy. In the TeV energy range bremsstrahlung photons become com- 
petitive with back-scattered laser photons in producing Be mesons. Although, the prospect 
of Be physics in 77 collisions are not very bright, observation of Be mesons does not appear 
completely unfeasible. 
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